Compound Poisson Statistics and Models of Clustering of Radiation Induced DNA 

Double Strand Breaks. 
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According to the experimental evidence damage induced by 
densely ionizing radiation in mammalian cells is distributed 
along the DNA molecule in the form of clusters. The most 
critical constituent of DNA damage are double-strand breaks 
(DSBs) which are formed when the breaks occur in both DNA 
strands and are directly opposite or separated by only a few 
base pairs. The paper discusses a model of clustered DSB 
formation viewed in terms of compound Poisson process along 
with the predictive assay of the formalism in application to 
experimental data. 
PACS numbers: 87.10. +e, 05.40. +j 



I. INTRODUCTION 

In living cells subjected to ionizing radiation many 
chemical reactions are induced leading to various bio- 
logical effects such as mutations, cell lethality or neo- 
plastic transformation The most important tar- 
get for radiation induced chemical transformation where 
these changes can be critical for cell survival is DNA dis- 
tributed within the cell's nucleus. Nuclear DNA is orga- 
nized in a hierarchy of structures which comprise the cel- 
lular chromatin. The latter is composed of DNA, histones 
and other structural proteins as well as polyamines. Or- 
ganization of DNA within the chromatin varies with the 
cell type and changes as the cell progresses through the 
cell cycle. Ionizing radiation produces variety of damage 
to DNA including base alterations and single- and double 
strand breaks (DSBs) in the sugar-phosphate backbone 
of the molecule Single strand breaks (SSBs) are 
efficiently repaired with high fidelity and probably con- 
tribute very little to the loss of function of living cells. 
On the other hand, DSBs are believed to be the critical 
lesions produced in chromosomes by radiation; interac- 
tion between DSBs can lead to cell killing, mutation or 
carcinogenesis. The purpose of theoretical modeling of 
radiation action [Q- is to describe qualitatively and 
quantitatively the results of radiobiological effects at the 
molecular, chromosomal and cellular level. The basic 
consideration in such an approach must be then descrip- 
tive analysis of breaks in DNA caused by charged particle 
tracks and by the chemical species produced. 
Production of DSBs in intracellular DNA can be stud- 
ied by use of the pulsed field gel electrophoresis (PFGE) 
B in which the gel electrophoresis is applied to elute 



high molecular weight DNA fragments from whole cel- 
lular DNA embedded in an organic gel (agarose). Two 
main approaches of this technique are usually applied. 
One is the measurement of the fraction of DNA leaving 
the well in PFGE, i.e. the amount of DNA smaller than 
a certain cutoff size defined by the electrophoretic condi- 
tions. This method has proven to be very sensitive, allow- 
ing reproducible measurements at relatively low doses. 
The second approach is to describe fragment-size distri- 
butions obtained after irradiation as a function of dose, 
taking advantage of the property of PFGE to separate 
DNA molecules based on how quickly they reorient in 
a switching (pulsed) electrical field. The major goal of 
the experiments is to quantify number of induced DSBs 
based on changes in the amount of DNA or the average 
fragment size in response to dose. In both cases data 
obtained are related to average number of DSBs. To an- 
alyze the data, the formalism describing random depolar- 
ization of polymers of finite size is usually adopted P|lC|] 
giving very well fits to experimental results with X-ray 
induced DNA fragmentation. In contrast to the findings 
for sparsely ionizing irradiation (X and 7 rays) charac- 
terized by low average energy deposition per unit track 
length (linear energy transfer, LETw 1 keV//zm), the 
densely ionizing (high LET) particle track is spatially lo- 
calized P,pJ|. In effect, multiplicity of ionizations within 
the track of heavy ions can produce clusters of DSBs 
on packed chromatin |ll|. The formation of clusters de- 
pends on chromatin geometry in the cell and radiation 
track structure. 

DSBs multiplicity and location on chromosomes may 
determine the distribution of DNA fragments detected 
in PFGE experiments. Modeling DNA fragment-size- 
distributions provides then a tool which allows to elu- 
cidate experimentally observed frequencies of fragments. 
Even without detailed information on the geometry of 
chromatin, models of radiation action on DNA can serve 
with some predictive information concerning measured 
DNA fragment-sizc-distribution. The purpose of the 
present paper is to discuss a model which can be used in 
analysis of DNA fragment-size- distribution after heavy 
ion irradiation. The background of the model is the Pois- 
son statistics of radiation events which lead to formation 
of clusters of DNA damage. The formation of breaks to 
DNA can be then described as the generalized or com- 
pound Poisson process for which the overall statistics of 
damage is an outcome of the random sum of random 
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variables (Section 2). Biologically relevant distributions 
are further derived and used (Section 3) in description of 
fragment size distribution in DNA after irradiation with 
heavy ions. Practical use of the formalism is discussed 
by fitting the distributions to experimental data. 



II. RANDOM SUMS OF RANDOM VARIABLES 
AND COMPOUND POISSON DISTRIBUTIONS 
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Consider [g_4 15| a sum Sn of N independent random 
variables X 



N 



Sn — ^ Xi 



(2.1) 



where N is a random variable with a probability gener- 
ating function g(s) 
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The above compound distribution is describing "clus- 
tered statistics" of events grouped in a number TV of 
clusters which itself has a distribution. As such, it is 
sometimes described in literature |]l6| as "mixture of dis- 
tributions" . Out of many interesting biological applica- 
tions of compound distributions P0|] , a special class 
constitute Poisson point processes which can be also an- 
alyzed in terms of random sums with Poisson distributed 
random events N. It can be shown that a mixture of 
Poisson distributions resulting from using any unimodal 
continuous function /(A) is a unimodal discrete distribu- 
tion. It is not so, however, in case of unimodal discrete 
mixing. In particular, mixtures of Poisson-Poisson or 
Poisson-binomial, known in literature as Neyman distri- 
butions can exhibit strongly multinomial character. 
By virt ue o f the above formalism and by using the for- 
mulae ( |2.5| ) , the generating function of the compound 
Poisson-Poisson distribution is: 
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g = exp(-A(l - /(«))) 



(2.7) 



and Xi are i.i.d. variables (independent and sampled 
from the same distribution) whose generating function 
/Wis 



where the random variables Xi are distributed according 
to a Poisson law 



f(s) = exp(-/i + fj,s) 



(2. 
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By use of the Bayes rule of conditional probabilities the 
probability that Sn takes value j can be then written as 



and the total Sn is a random variable with a compound 
Poisson-Poisson (Neyman type A) distribution: 

^ (Nu) x e~ Nfl \ N e~ x 
P(S N = x)=P(x;»,\) = ^ [ -^Ll ^JL_ (2 .9) 
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P(S N = j) = hj = V P(S N = j\N = n)P(N = n) (2.4) for which the mean and variance are given by 



n=0 

For fixed value of n and by using the statistical inde- 
pendence of Xi's, the sum Sn has a probability gener- 
ating function F(s) being a direct product of /(s), i.e. 
F(s) = f( s ) n = Yjf=Q F j^ from which it follows that 
P(S N = j\N = n) = Fj. The formula (|J) leads then 
to the compound probability generating function of Sn 
given by 



E[x] = /iA, Var[x] = A/i(l + /x) 



(2.10) 



h(s) = hjsi = 



3=0 



oo 

= ^ S „/( S )"e 9 {/( S )} 

n=0 



(2.5) 



Conditional expectations rules can be used to determine 
moments of a random sum. Given E[N] = v, E\X,\\ = /x, 
T^ar[iV] = r 2 and yar[ATi] = a 2 , the first and the second 
moment of the random sum Sn are 



The resulting distribution can be interpreted as a mix- 
ture of Poisson distribution with parameter N/i where 
N (number of clusters) is itself Poisson dist ribu ted with 
parameter A. Figures 1,2 present function ( |2.9| ) for two 
various sets of parameters A, /i. 

The compound Poisson distribution (CPD) has a wide 
application in ecology, nuclear chain reactions and que- 
ing theory PJl9|-[2l]] . It is sometimes known as the distri- 
bution of a "branching process" and as such has been 
also used to describe radiobiological effects of densely 
ionizing radiation in cells [ |17|j22| - |24]| . When a single 
heavy ion crosses a cell nucleus, it may produce DNA 
strand breaks and chromatin scissions wherever the ion- 
izing track structure overlaps chromatin structure. The 
multiple yield of such lesions depends on the radial distri- 
bution of deposited energy and on the microdistribution 
of DNA in the cell nucleus. The latter and the geome- 
try of DNA coiling in the cell nucleus determine number 
of crossings, the "primary" incidents leading to DSBs 
production. By assuming for a given cell line, a "typi- 
cal" average number n of possible crossings per particle 
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traversal, the distribution of the number of chromatin 
breaks i can be modelled by a binomial law: 



P(i\n) 



(2-11) 



where p is a probability that a chromatin break occurs 
at each particle crossing (and q is the probability that 
it does not). The overall probability that i lesions will 
be observed after m independent particles traversed the 
nucleus is given by M 
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which is a compound Neyman type B distribution ob- 
tained as a random Poisson sum of binomially distributed 
i.i.d variables. In the above presentation the average 
number of particles crossing the cell nucleus A is pro- 
portional to the absorbed energy (dose) and given by a 
product A = <tF of particle fluence F and nuclear cross 
section a. 



P N (0\a,F,n) = 

m—l 



(nmy.q nm (crF) m e- aF _ 
(nm)\m\ 

= exp[-aF(l - q n )} (2.14) 



which differs by a factor (1 — q n ) in the exponent from 
the surviving fraction for a Poisson distribution: 



P P (0\a,F,n) = exp[-cr^] = cxp[-E[i] 



(2.15) 
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FIG. 1. Simulated probability density function for the Ney- 
man-type A distribution (^1]) with \ — 6, /i — 100 for 
N — 10000 points. Note the finite value at x = corre- 
sponding to P(0; n, A). 

Aggregation of observed cellular damage potentially 
leads to the phenomenon of "overdispersion" - that 
is, the variance of the aggregate may be larger than 
Poisson variance yielding "relative variance" Var re i = 
Var[SN]/ E[Sn] larger than 1. Assuming thus the Pois- 
son statistics of radiative events, for any distribution of 
lesions per particle traversal, the condition for overdis- 
persion can be easily rephrased in terms of (2.6) 



Var[Xi]/E[Xi] + E[X t \ < 1 



(2.13) 



If no repair process is involved in diminishing number of 
initially produced lesions, the sur vivin g fraction of cells 
can be estimated from formula eq.(2.12) as a zero class of 
the initial distribution, i.e. the proportion of cells with 
no breaks 



FIG. 2. Simulated probability density function for the Ney- 
man-type A distribution ( [2.9| ) with A = 100, fj, = 6 for 
N = 10000 points 



III. DNA FRAGMENTS DISTRIBUTION 
GENERATED BY IRRADIATION: STATISTICAL 
MODEL. 

DNA double stranded molecules in a size range from 
a few tenths of kilobase pairs to several megabase pairs 
can be evaluated by the PFGE technique. Randomly dis- 
tributed DSBs are detected as smears of DNA fragments. 
The DNA mobility mass distribution may be transformed 
into a fragment length distribution using a calibration 
curve. It is obtained by relating migration distance of 
DNA within the gel to molecular length with the aid of 
size markers loaded on the same gel p5[ . To interpret 
the experimental material one needs to relate percentage 
of fragments in defined size ranges to number of induced 
DSBs. For that purpose several models have been de- 
rived, mainly based on the description of random depo- 
larization of polymers of finite size [p| JlO|j26]] . Although 
the models give satisfactory prediction of size-frequency 
distribution of fragments after sparsely ionizing radiation 
[i.e for X-rays and 7), they generally fail to describe the 
data after densely ionizing radiation 13 2j|. The exper- 
iments with heavy ions demonstrate that after exposure 
to densely ionizing particles gives rise to substantially 
overdispersed distribution of DNA fragments which indi- 
cates the occurrence of clusters of damage. The following 
analysis presents a model which takes into account forma- 
tion of aggregates of lesions after heavy ion irradiation. 
Fragment distribution in PFGE studies is measured by 
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use of fluorescence technique or radioactive labeling with 
the result being the intensity distribution. The generated 
signal is proportional to the relative intensity distribution 
of DNA fragments and can be expressed as 



with 



D(x) 



I(x) = xD(x) 



3=0 



(3.1) 



(3.2) 



where D(x\j) stands for the density of fragments of 
length x provided j DSBs occur on the chromosome of 
size S. Frequency distribution of the number of DSBs is 
assumed here in the form of CPD (|j]) with parameters fx 
and A representing average number of breaks produced by 
a single particle traversal and average number of particle 
traversals, respectively. The "broken-stick" distribution 
p7| , ^6| for j breaks on a chromosome of size S yields a 
density of fragments of size x: 



D(x\j) = 8(x-S) + 2j±(l-l; 



Did 



s 



3-1 



(3.3) 



where the first two terms describe contributions from the 
edge fragments of the chromosome and the third term de- 
scribes contribution from the internal fragments of length 
x < S. The first term applies to the situation when 
j = 0; the edge contribution can be understood by ob- 
serving that the first and the j + 1 fragment have the 
same probability of being size x. Direct summation in 
formula (|3.2|) leads to 



D N (x) = exp(-A(l - e-^))5{x - S) + 
A(e-"#-l)) + 



2A/i x 
— exp(-„- 



-e- A (l - |)^(1 + Ae""S ) exp(- M | + Ae""*) (3.4) 



for Neyman distribution of number of breaks j and to 

(3.5) 
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for a Poisson distribution with parameter A. Integration 
of I(x) (eq.(3.1)) from to some average (marker) size 
X* and division by S yields the relative fraction of DNA 
content. For A >> 1 and fi << 1, the Neyman- type 
A distribution converges to a simple Poisson. In such a 
case, simplified expression (3.4) leads to results known 
in literature as "Blocher formalism" [p| Jir^j2r| which de- 
scribes well the DNA content in probes irradiated with 
X- and 7-rays. 
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FIG. 3. Distribution of DNA content (integrated eq.(fu])) 
as a function of the dose and fragment size for 
S — 245Mbp, (i — 5. The fragments length is in Mbp units. 

Figure 3 presents predicted dose-response curves for 
the model. The amount of DNA content is shown in 
function of dose and fragment size. In calculation, the 
parameter S — 245 mega base pairs has been used which 
is the mean chromosome size for Chinese hamster cells, 
the cell line for which experimental data are displayed in 
Figure 4. 

The increase in multiplicity of DSBs produced per one 
traversal of a particle leads to pronounced increase in 
production of shorter fragments which is illustrated in 
the shift of the peak intensity towards smaller x values. 
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FIG. 4. Fraction of DNA content observed experimentally 
within the range of sizes 0.1-1.0 Mbp. Data show higher prob- 
ability of producing short fragments after irradiation with 
particles than for sparsely ionizing radiation at comparative 
dose. Lines represent the best fit to eq.(3.1) by use of D^(x) 
function for heavy ions (Au: A = 3 x 10~ 3 ,/i = 6x 10 2 ; C: 
A = 6 x 10~ 3 ,^i = 6 x 10 2 ) and D P {x) for X-rays (A = 0.85). 
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IV. SPATIAL CLUSTERING OF BREAKS AND 
NON-POISSON STATISTICS. 

Clustering of breakage events can be viewed as the pro- 
cess leading to non-exponential "spacing" between sub- 
sequent events, similar to the standard analysis of level 
repulsion in spectra of polyatomic molecules and com- 
plex nuclei. For a random sequence, the probability that 
a DSB will be in the infinitesimal interval 



(X + x, X + x + dx) 



(4.1) 



proportional to dx is independent of whether or not there 
is a break at X. This result can be easily changed by 
using the concept of breaks "repulsion'. Given a break 
at X, let P(x)dx be the probability that the next break 
(x > 0) be found in the interval (X + x, X + x + dx). We 
then have for the nearest-neighbour spacing distribution 
of breaks the following formula: 



P(x)dx = Prob(l G dx\0 G x)Prob(0 G x) 



(4.2) 



where Prob(n G dx\m G x) is the conditional probability 
that the infinitesimal interval of length dx contains n 
breaks wheras that of length x contains m of those. The 
first term on the right-hand side of the above equation 
is dx times a function of x which we denote by r(x), 
depending explicitly on the choices 1 and of the discrete 
variables n and m. The second term is given by the 
probability that the spacing is larger than x: 



(4.3) 



P(y)dy 
Accordingly, one obtains 

P(x) = r(x) / P{y)dy, 

J X 

whose solution can be easily found to be 

r 

P{x) = Cr(x) exp(— / r(y)dy) 



(4.4) 



(4.5) 



where C is a constant. The Poisson law, which reflects 
lack of correlation between breaks, follows if one takes 
r(x) — A, where A -1 is the mean spacing between DSBs. 
If choosing on the other hand 



r(x) = \x 



A-l 



(4.6) 



i.e. by assuming clustering of points (DSBs) along a line, 
one ends up with the Weibull density. The constants C 
and A can then be determined from appropriate condi- 
tions, e.g. 



Jp(x)dx = l, 



(4.7) 



xP{x)dx = A" 1 

One then finds that 

P(x) = \e~ Xx 
for the Poisson distribution and 

P{x) = Xx*- 1 cxp(-x x ) 



(4.8) 



(4.9) 



(4.10) 



for the Weibull analogue. Note that the above de nsit y 
can be derived as a generalization of the law eq.( |4.9| ): 
the Weibull density can be obtained as the density of 
random variable y = x x / x with x being an exponential 
random variable. For A > 1, the Weibull distribution is 
unimodal with a maximum at point x m = (1 — A _1 ) A 
In this one easily recognizes for A = 2 the spacing dis- 
tribution of the Wigner law. The latter displays "re- 
pulsion" of spacing, since P(0) = 0, in contrast to the 
Poisson case which gives maximum at x = 0. Fractional 
exponent A < 1 describes, on the other hand, enhanced 
frequency of short spacings which, in fact, matches bet- 
ter experimental data for heavy ions(c/. Figure 4). The 
above analysis brings also similarities with random walks 
p9| , p0[ where symmetry breaking transition manifests it- 
self as a change in the spectral spacing statistics of decay 
rates. In such cases, the statistics of events of interest de- 
viates, as a counting process, from the regularity of Pois- 
son process, for which the subsequent event arrivals are 
spaced with a constant mean A . The clustered statis- 
tics of breakage can be thus viewed as a (fractal) random 
walk or a cumulative distribution of a random sum of 
random variables eq.(2.1). The problem of characteriz- 
ing the limit distribution for such cases with underlying 
"broad" distributions g(x) of Xi has been studied ex- 
tensively in mathematical literature [ [32"| and has been 
solved with classification of the possible limit distribu- 
tions provided that requirement of "stability" is fulfilled 
under convolution. Following the definition, the distri- 
bution g(x) is stable, if for any N there exist constants 
cjy > and d/v such that Sn has the same density as the 
variable y = c^Xi + d^. The stability condition can be 
rephrased in terms of the canonical representation given 
by a form of the characteristic function (i.e. the Fourier 
transform g(k)) of stable distributions I33U32] 



\ng(k) = ijk - C|fc| A [l - iuj(3sign(k)} 



(4.11) 



where 7 is real, < A < 2, 10 is real and |w| < tan(irX/2)\. 
The cases relevant for biological modelling are covered 
by 1 < A < 2 (stable distributions have no variance if 
A < 2 and no mean if A < 1). In particular, positiv- 
ity of steps in the random walk modelled by eq.(2.1) al- 
lows for g(k) = exp[~C\k\ x ] which gives asymptotically 



-A-l 



Probability distribution that x > z sat- 



and 



isfies then f(z) ~ z for z — > 00. The resulting distri- 
bution is "self-similar" in the sense that rescaling z to 
Az and f(z) to A~ x f(z) does not change the power law 



5 



distribution. In other words, the number of realizations 
larger than Az is times the number of realizations 
larger than z. The power-law probability distribution 
function describes then the same proportion of shorter 
and larger fragments whatever size is discussed within 
the power law range. For A = 1/2, C = 1, u> = 1 the form 
of Levy-Smirnov law is recovered 
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g(x) = (2n)- 1 / 2 x- 



-3/2 e -A 



(4.12) 



The probability density eq.( 1.12] ) has a simple interpre- 
tation as the limiting law of return times to the origin for 
a one-dimensional symmetrical random walk and as such 
has been also used to describe the fragment size distribu- 
tion of a one dimensional polymer pl|j3^] . In the prob- 
lems related to polymer fragmentation induced by irradi- 
ation, the approach based on a random walk with fluctu- 
ating number of steps (or, equivalcntly, on a point proce- 
ses model with a clustered statistics of waiting times) is 
a legitimate one as it can comprise the natural random- 
ness of primary events (i.e. particle hits of biological 
target) and secondary induction of multiple (clustered) 
lesions. Further investigations in this field should lead to 
better understanding of possible emergence of power-law 
distributions of larger fragments on kbp and Mbp scales. 



V. CONCLUSIONS 

An existing substantial evidence demonstrates that ex- 
posure to densely ionizing charged particles gives rise to 
overdispersed distribution of chromatin breaks and DNA 
fragments which is indicative of clustered damage occur- 
ing in irradiated cells. The clustering process can be 
expressed for any particular class of events such as ion- 
izations or radical species formation and is a consequence 
of energy localization in the radiation track. Chromoso- 
mal aberrations expressed in irradiated cells are formed 
in process of misrejoining of fragments which result from 
production of double-strand breaks in DNA. The loca- 
tion of double-strand breaks along chromosomes deter- 
mines DNA fragment-size distribution which can be ob- 
served experimentally. The task of stochastic modeling 
is then to relate parameters of such distributions to rele- 
vant quantities describing number of induced DSBs. Ap- 
plication of the formalism of clustered breakage offers 
thus a tool in evaluation of the radiation respone of DNA 
fragment-size distribution and assessment of radiation in- 
duced biological damage. 



Acknowledgements . 

E.G-N acknowledges partial support by KBN grant 2P03 
98 14 and by KBN-British Council collaboration grant 
C51. 



[1] E.L. Alpen, Radiation Biophysics, (Academic Press, San 

Diego, 1998). 
[2] G.Kraft, Nucl. Science Appl. 1 (1987) 1. 
[3] C. Von Sonntag, The Chemical Basis of Radiation Biol- 
ogy, (Taylor and Francis, London, 1987); J.F. Ward, Int. 

J. Radiat. Biol. 66 (1994) 427. 
[4] C.A. Tobias, E. Goodwin and E. Blakely, in Quantitative 

Mathematical Models in Radiation Biology, J. Kiefer, ed., 

Springer Verlag, Berlin 1988, p. 135. 
[5] P.J. Hahnfeldt, R.K. Sachs and L.R. Hlatky, J. Math. 

Biol. 30 (1992) 493. 
[6] A.Chatterjee and W. Holley, Int. J. Quant. Chem. 391 

(1991) 709. 

[7] R. Sachs, D.J. Brenner, P.J. Hahnfeldt and R. Hlatky, 

Int. J. Radiat. Biol. 74 (1998) 185. 
[8] G. Iliakis, D. Blocher, L. Metzger and G. Pantelias, Int. 

J. Radiat. Biol. 59 (1991) 927. 
[9] E.W. Montroll and R. Simha, J. Chem. Phys. 8 (1940) 

721. 

[10] D. Blocher, In. J. Radiat. Biol. 57 (1990) 7. 
[11] M. Kramer and G. Kraft, Rad. Env. Biophysics 33 (1994) 
91. 

[12] G. Taucher-Scholz and G. Kraft, Rad. Res. 151 (1999) 
595 

[13] M. Lobrich, P. Cooper and B. Rydberg, Int. J. Ra- 
diat. Biol. 70 (1996) 493; H.C. Newman, K.M. Prise, 
M. Folkard and B.D. Michael, ibid 71 (1997) 347; E. 
Hoglund, E. Blomquist, J. Carlsson and B. Sternlow, ibid 
76 (2000) 539. 

[14] N.G. Van Kampen, Stochastic Processes in Physics and 
Chemistry, (North Holland, Amsterdam, 1981). 

[15] A. Papoulis, Probability, Random Variables and Stochas- 
tic Processes, (McGraw-Hill, Tokyo, 1981). 

[16] M. Kendall and A. Stuart, The Advanced Theory of 
Statistics, Charles Griffin & Co., London, 1977. 

[17] N. Goel and N. Richter-Dyn, Stochastic Processes in Bi- 
ology, (Academic Press, New York, 1974). 

[18] T. Maruyama, Mathematical Modeling in Genetics, 
(Springer Verlag, Berlin, 1981). 

[19] AT. Bharucha-Reid, Elements of the Theory of Markov 
Processes and Their Applications, (Dover Publications, 
New York, 1988). 

[20] S. Karlin and H. Taylor, First Course in Stochastic Pro- 
cesses, (Academic Press, New York, 1976). 

[21] J. Neyman, Am. Math. Stat. 10 (1939) 35. 

[22] N. Albright, Radiat. Res. 118 (1989) 1. 

[23] E. Gudowska-Nowak, S. Ritter, G. Taucher-Scholz and 
G. Kraft, Acta Phys. Pol. 31B (2000) 1109. 

[24] E. Nasonova, E. Gudowska-Nowak, S. Ritter and G. 
Kraft, Int. J. Radiat. Biol., (2000), in press. 

[25] J. Heilmann, G. Taucher-Scholz and G.Kraft, Int. J. Ra- 
diat. Biol., 68 (1995) 153; G. Taucher-Scholz and G. 
Kraft, Rad. Res. 151 (1999) 595 



G 



[26] T. Radivoyevitch and B. Cedervall, Electrophoresis 17 
(1996) 1087. 

[27] P.J. Flory Statistical Mechanics of Chain Molecules (In- 

terscience, New York, 1969). 
[28] G. Van den Engh, R. Sachs and B. Trask, Science 257 

(1986) 1410. 

[29] P. Alpatov and L.E. Reichl, Phys. Rev. E 52 (1995) 4516. 
[30] D.R. Nelson and N.M. Shnerb, Phys. Rev. E. 58 (1998) 
1384. 

[31] G.H. Weiss and R.J. Rubin, Adv. Chem. Phys. 52 (1983) 
363. 

[32] V.M. Zolotarev One-dimensional Stable Distributions, 
(American Mathematical Society, Providence, 1986) 

[33] B.V. Gnedenko and A.N. Kolmogorov Limit Distribu- 
tions for Sums of Independent Random Variables, (Ad- 
dison Wisley MA 1954) 

[34] A.L. Ponomarev and R.K. Sachs, Bioinformatics 15 
(1999) 957. 



7 



